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Classical correlations of ground states typically decay exponentially and polynomially, respectively 
for gapped and gapless short-ranged quantum spin systems. In such systems, entanglement decays 
exponentially even at the quantum critical points. However, quantum discord, an information- 
theoretic quantum correlation measure, survives long lattice distances. We investigate the effects of 
quenched disorder on quantum correlation lengths of quenched averaged entanglement and quantum 
discord, in the anisotropic XV and XYZ spin glass and random field chains. We find that there is 
virtually neither reduction nor enhancement in entanglement length while quantum discord length 
increases significantly with the introduction of the quenched disorder. 


I. INTRODUCTION 

In the last twenty years or so, it has undoubtedly 
been established that entanglement [1] plays a significant 
role in efficient quantum communication protocols [2], 
which include quantum dense coding [3], quantum tele¬ 
portation [4] and quantum key distribution [5], and in 
quantum computational tasks like the one-way quantum 
computer [6]. Generation of entanglement in compos¬ 
ite systems, requires interaction between its subsystems 
and hence strongly interacting systems, such as quantum 
spin chains, form natural resources of entanglement. In 
the last decade, it has also been realized that entangle¬ 
ment can be used to detect low-temperature phenomena 
in many-body systems [7] . Spin models have further been 
identified as ideal setups for transferring quantum states 
[8] as well as for realizing quantum gates, essential parts 
of quantum circuits and quantum simulators [6, 9]. More¬ 
over, spin models can be prepared in laboratories via e.g. 
ultracold atoms in optical lattices [10], trapped ions [11], 
and nuclear magnetic resonance systems [12] in a con¬ 
trolled way, and also can be found in materials [13, 14]. 

Distribution of entanglement among various parts of 
a many-body system is important for transmitting in¬ 
formation from one party to another within communica¬ 
tion computing networks [6, 8, 15]. The efforts to under¬ 
stand the sustainability of correlations in many body sys¬ 
tems as a function of lattice distance happens to be cru¬ 
cial, since the low-temperature scaling of correlations are 
expected to remain universal, irrespective of the micro¬ 
scopic details of various materials governed by the same 
class of underlying models. In case of gapped systems 
with short-ranged interactions, using the Lieb-Robinson 
bounds [16], it can be shown that two-site classical corre¬ 
lation functions, {OiOj) — {Oi){Oj), decay exponentially 
with increasing lattice distance, known as clustering of 
correlations [17, 18] except when the systems are critical, 
when they decay polynomially, with Oi and Oj being 
the observables at sites i and j respectively. However, 
it has been reported that in many cases, quantum corre¬ 
lations, as qualified by entanglement, decays much more 
rapidly in comparison to classical correlators, both near 
and far from the quantum critical points. For example. 


it was shown that in the quantum XY model with trans¬ 
verse magnetic field, entanglement can survive only up 
to the next nearest neighbor [19], while in the Heisen¬ 
berg model, it becomes vanishingly small already after 
the nearest neighbor [20]. 

In recent times, information-theoretic quantum corre¬ 
lation measures, like quantum discord [21] and quan¬ 
tum work deficit [22], were introduced which promise 
to estimate quantum correlations beyond entanglement 
[23]. For example, it has been shown that there exist in¬ 
stances in which quantum discord is capable of detecting 
quantum phase transition while pairwise entanglement 
fails, especially at finite temperature [24, 25]. Hence, 
it is natural to investigate the scaling of these two-site 
information-theoretic quantum correlation measures in 
these systems, with increasing distance between the sites. 
In this respect, it is already known that unlike entangle¬ 
ment, quantum discord can survive over a longer distance 
[26, 27] in the ordered XY spin chain. 

Establishing finite correlations between distant parties 
is undoubtedly important to implement quantum infor¬ 
mation processing tasks in many body systems. In this 
paper, we ask the following question: Is it possible to en¬ 
hance quantum correlation lengths, namely entanglement 
and quantum discord lengths, significantly, by introduc¬ 
ing defects in quantum spin systems? We find that the 
answer is in the affirmative. Disorder occurs unavoidably 
in real materials and can now also be engineered artifi¬ 
cially in, e.g. cold atom experiments [28]. Moreover, it 
was shown that there exists some models, in which disor¬ 
der plays a constructive role by enhancing physical prop¬ 
erties like magnetization, classical correlators, and quan¬ 
tum correlations [14, 29, 30]. Such phenomena, known as 
“order-from-disorder” or “disorder-induced order”, run 
certainly in contrary to the naive belief that impurity in 
the systems can have only a debilitating effect. 

In this work, we investigate quantum correlation 
lengths for the ground states of one-dimensional 
quenched disordered anisotropic XY and XYZ spin 
chains. Our results show that although the entangle¬ 
ment lengths cannot be improved significantly in these 
quenched disordered systems, the lengths of information- 
theoretic quantum correlation measures, can substan¬ 
tially be enhanced in such systems. Such disorder- 
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induced advantages can be seen in quantum XY spin 
glass and random field XY models, where the investi¬ 
gations are carried out with the help of Jordan-Wigner 
transformations, and in the XYZ spin glass systems, 
which are handled by using density matrix renormaliza¬ 
tion group (DMRG) techniques. 

The rest of the paper is arranged as follows. Section II 
introduces the models under study, and discusses the 
methods involved in solving them. Section III reviews 
various quantum correlation measures used in this work 
and introduces quantum correlation lengths. The results 
for the disordered XY and XYZ models are discussed in 
Secs. IV and V respectively. Finally, we summarize our 
results in Sec. VI. 


II. MODELS AND METHODOLOGY 

In this section, we briefly describe the models that we 
use in this paper. Among the models described, note 
that the ordered XY spin model is exactly solvable and 
the analytic technique can be used to handle the corre¬ 
sponding disordered model, while the XYZ spin model 
with and without disorder cannot be solved analytically. 
We also present a brief decription of quenched averaging. 


A. Quenched disorder and averaging 

The disorder in the system parameters are taken to be 
“quenched”. That is, we assume that the time scale over 
which the system dynamics of interest takes place is much 
smaller compared to the time scale over which there is 
a change in the particular set of parameters governing 
the disorder in the system. In order to calculate the 
quenched averaged value of a physical quantity, we need 
to perform the averaging over the probability distribution 
of several realizations, each of which corresponds to a 
fixed configuration of the system, after calculating the 
value of the physical quantity for the fixed configurations. 


B. Quantum XY spin chain: Ordered and 
Disordered Models 


The general Hamiltonian for the quantum XY spin 
chain with nearest-neighbor interactions in an external 
magnetic field is given by 


^ = ""[Zl 7 ((1 + 7)«+i + (1 - 7 ) 0 - 1 'crf+i) 



where kJ^ are the coupling constants, nhi is the magnetic 
field strength at the site, and 7 is the anisotropy con¬ 
stant. The constant k has the units of energy, while Ji^hi^ 


and 7 are dimensionless. Here, cr^, for j = x,y^z, corre¬ 
spond to the Pauli spin matrices. Moreover, we assume 
the periodic boundary condition, i.e., ctat+i = ai. 

Case 1: Quantum XY spin glass:: In this case, 
the coupling strengths Ji are randomly chosen from in¬ 
dependently and identically distributed (i.i.d.) Gaussian 
distribution with mean (J) and unit standard deviation. 
However, the field is kept uniform throughout the lattice, 
i.e., hi = h for i = 1, • • • , N. 

Case 2: Random field quantum XY spin chain:: 
The model is now with uniform coupling, i.e., 3% = J but 
the hi are i.i.d. Gaussian random variables with mean 
{h) and unit standard deviation. 

Case 3: The ordered quantum XY spin chain:: 
The model is with site-independent coupling constants 
as well as the magnetic field strengths, i.e., Ji = J and 
hi = h for i = 1, • • • ,7V. 

The ordered quantum XY spin chain is exactly solv¬ 
able via successive applications of the Jordan-Wigner, 
the Fourier, and the Bogoliubov transformations [31-33]. 
The corresponding disordered systems can also be han¬ 
dled upto relatively large system sizes by using the same 
transformations. The two-site reduced density matrix, 
Pij, of the ground state can be easily constructed from 
the one and two-point correlation functions as 



0/) -kTO|(/0cr^) 


OL—x,y,z 


( 2 ) 


Note that mf = = 0 and all the off-diagonal correla¬ 

tions vanish. 


C. Quantum XYZ spin glass 


Here, our interest lies in the quenched averaged corre¬ 
lation lengths in the XYZ spin chain. The Hamiltonian 
of the quantum XYZ spin glass model is given by 


N-l 


H = K 


[ Z [j [(1 + 7)«+i + (1 - + 


2=1 


A , , 

<^i+i 




where kA is the nearest-neighbor coupling strength for 
the zz- interaction, which is independent of the site. The 
rest of the parameters are the same as discussed above 
in the context of the XY spin glass. The corresponding 
ordered Hamiltonian can be obtained from Eq. (3) by 
simply setting Ji = J toi i = 1,..., N. 

Unlike the quantum XY disordered chain, for which 
the ground state of a considerably large number of spins 
can be obtained by using the Jordan-Wigner transforma¬ 
tion, one needs to resort to numerical techniques for the 
XYZ model with random coupling strengths. In order 
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to investigate the ground state for the system charac¬ 
terized by the Hamiltonian in Eq. (3), we employ the 
well-established numerical technique called the DMRG 
method [34]. After performing the standard infinite size 
DMRG, several finite size DMRG sweeps are also carried 
out in order to increase the accuracy of the calculations 
for the inhomogeneous chain. 


III. QUANTUM CORRELATION MEASURES 

The presence of quantum correlations between subsys¬ 
tems of a composite system helps in realizing many quan¬ 
tum information protocols. In order to explore these 
protocols, it is necessary to quantify the quantum cor¬ 
relations involved. In this work, we have mainly used 
two quantum correlations measures, namely, concurrence 
and quantum discord. They belong to two different 
paradigms of quantum correlation - while the first cor¬ 
responds to the entanglement-separability paradigm, the 
other is in the information-theoretic one. In the follow¬ 
ing subsections, we will briefly introduce both the quan¬ 
tum correlation measures considered here. We will then 
provide a short introduction of the concept of quantum 
correlation lemgth. 


A. Concurrence 

Concurrence [35] quantifies the amount of entangle¬ 
ment present in an arbitrary two-qubit state. Given a 
two-qubit density matrix, pab-, the concurrence is de¬ 
fined as 

C{pab) = max{ 0 , Ai - A2 - A3 - A4}, (4) 

where A^’s are the eigenvalues of the Hermitian matrix 
R = \/i/]o p and satisfy the order Ai > A 2 > A 3 > A 4 . 
Here p = (ay ^ ay) p* (ay ^ ay) with p* being the complex 
conjugate of p in the computational basis. 


B. Quantum discord 

In classical information theory, the amount of igno¬ 
rance about a probability distribution, is quanti¬ 
fied by the Shannon entropy, defined as H({pi}) = 
— Pi log2Pi- The mutual information between two 
classical random variables i and j, having the marginal 
distributions {pi} and {pj}-, can be defined in two equiv¬ 
alent ways as 

AiPij}) = H({p^}) +H({pj}) - H({p,j}) 

where {pij} and {Pi\j} correspond to the joint probability 
distribution of the variables i and j, and the conditional 


probability distribution respectively. In case of quantum 
systems, the quantity 

I(pab) = S(pa) + S(pb) - S(pab), ( 6 ) 

for a two-party quantum system pab, can be argued 
to quantify the total correlation present in the system 
where S(a) = —tr[(T log 2 cr] and p^, i = A^B are the lo¬ 
cal density matrices of pab [36-38]. This quantity can 
be interpreted as the quantized version of the first ex¬ 
pression of the classical mutual information in Eq. (5). 
To quantize the other expression, we consider that a 
measurement is performed on one party, say B, using 
a complete set of rank-one projectors, {Bi}, satisfying 
the relations BiBj = StjBi and J2iBi = Ib- The post¬ 
measurement ensemble is given by {pi , p\^} with p\^ = 

(lA®Bi) p (lA®Bi) andp* = tT:((lA®Bi) pab (lA®Bi)), 

where I a is the identity operator on the Hilbert space of 
the system with observer A. The corresponding quantum 
conditional entropy is given by 

<S'(pa|b) = min^PiS'(p^li), (7) 

^ i 

where pA\i = ^t^bPab- quantum version of the sec¬ 
ond expression in Eq. (5) of classical mutual information 
then reads 

J(pab) = S(pa) - S(paib), ( 8 ) 

which turns out to be inequivalent to the expression in 
Eq.(6), and can be argued as a measure of classical corre¬ 
lation of the state pab- Quantum discord [21] is defined 
as the difference between these two inequivalent quanti¬ 
ties, and is given by 

T^(pab) =I(pab) — J(pab)- ( 9 ) 

In all the cases considered in this paper, the bipartite 
states are X states [39] with for which 

an analytical form of quantum discord is available [40- 
43]. We have also checked the claim numerically. When 
the measurement is performed by the first party, i.e. on 
the spin of a two-party state pij , the classical correla¬ 
tion is given by J(pi,j) = H 2 - H 2 (^), with 

p = ^J(ml)‘^ + (TfJ)^ and H 2 (x) = -x\og 2 X - (1 - 
x) log 2 (l — x) being the binary entropy. 

C. Quantum correlation length 

As mentioned earlier, our primary aim is to investi¬ 
gate quantum correlation length in quenched disordered 
XY and XYZ spin models. For arbitrary ordered spin 
systems, let us first define the entanglement length, in 
particular, concurrence length. If the concurrence, Cij 
between the and spins behaves as 

li-il 

Cij ~ e , 


( 10 ) 
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we refer to as the concurrence length. For quenched 
disordered systems, the concurrence length for the 
quenched averaged two-site concurrences, if the latter be¬ 
haves as in Eq. (10), will be denoted by ^{c>- In this pa¬ 
per, entanglement length will mean concurrence length. 

Similar to the concurrence length, one can define dis¬ 
cord length, for the ordered models, if the quantum 
discord, Dij, between the and the sites behaves 
as 

Dij ^ a + b e (11) 

where a, b are parameters that are determined by the 
model under consideration. And similarly as above, the 
discord length for quenched disordered spin models can 
be defined, and is denoted by ^{d)- A note on the units 
used here and henceforth in the paper. The entangle¬ 
ment and discord length are in the units of lattice dis¬ 
tance. The parameters a and b in Eq. (11) and in similar 
equations are in bits, just as for quantum discord. We 
have kept silent a multiplicative constant on the right 
hand side of Eq. (10) that has the unit of ebits, just as 
the concurrences. 

Note that this is different from the concept of localiz- 
able entanglement, the amount of entanglement that can 
be concentrated between two parties through measure¬ 
ments performed on the rest of the parties [44] (see also 
[45-47]). 

IV. ENTANGLEMENT AND DISCORD 
LENGTHS IN ANISOTROPIC XY MODEL 

The ordered quantum XY model is exactly solvable 
for hnite as well as infinite spin chains [31-33]. In the 
thermodynamic limit of the system, closed form expres¬ 
sions for the magnetization [32] and the two-point corre¬ 
lation functions [33] are known, which are necessary to 
calculate the two-body density matrices between two ar¬ 
bitrary sites, and eventually to compute quantities like 
the concurrence and quantum discord. However, this is 
not the case for the disordered systems and one is ulti¬ 
mately restricted to study finite size systems, but tech¬ 
niques similar to those for the infinite system can help 
to compute quantum correlations for relatively large sys¬ 
tems [31, 48, 49]. 

In order to make a comparison between the disordered 
and ordered systems, we fix the means of the distribu¬ 
tions of the disordered parameters in the disordered sys¬ 
tems to be identical to the corresponding parameters of 
the ordered system. In particular, the quenched averaged 
physical quantity, Qavi{a), (b),...), corresponds to the 
disordered system with disordered parameters a,b,... 
having means (a) ,{b),... respectively. The correspond¬ 
ing physical quantity for the ordered system is then 
Q{{a),{b ),..where the values of the system parame¬ 
ters a,b,... are kept constant at (a), (&),... respectively. 

Note that in the case of the ordered XY model, the 
ground state is multipartite entangled except at J/h = 
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FIG. 1. (Color online.) Entanglement length in qnantnm XY 
spin glass vs. the ordered XY model. In each panel, the dot¬ 
ted curve denotes the quenched averaged concurrence, (Cij), 
between the sites i and j, plotted against the lattice distance 
r = \i — j\, for the quantum XY spin glass for a certain value 
of {J)/h. The corresponding curves for the ordered systems 
are also drawn in each panel. The vertical axis denote the 
concurrences, while the horizontal ones denote the lattice dis¬ 
tances. Here, N = 50 and 7 = 0.5. For the disordered case, 
the number of realizations of the random coupling is taken to 
be 10^, for the quenched averaging. The lines are exponential 
fits and the obtained data at integer values of r. The vertical 
axis are measured in ebits, while the horizontal axis are in 
lattice length unit. 

l/x/l" — 7 ^ which is known as the factorization point 
[50, 51]. At this point, the ground state is doubly degen¬ 
erate and both the degenerate states are factorized as a 
tensor product of quantum states corresponding to all the 
individual spins. Since we have taken the ground state 
to be a pure symmetric state by taking equal superposi¬ 
tion of both the degenerate ground states, the two spin 
entanglement here vanishes at this point, while quantum 
discord may have a non-zero value. 


A. Anisotropic quantum XY spin glass: 
Entanglement length vs. Discord length 

Let us first consider the quantum XY spin glass with 
N sites. All considerations are for the pure symmetric 
ground state. 


1. Entanglement length 

We will now investigate the effects of impurities on en¬ 
tanglement length in the XY spin model. To this end, we 
compare the two-site concurrence between the sites i and 
j of the ordered system, described by H{J,h), with the 
quenched two-site concurrence between the same sites for 
the XY spin glass system with {J)/h = J/h, where {J) is 
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FIG. 2. (Color online.) All considerations here are the same 
as in Fig. 2, except that quantum discord is considered instead 
of concurrence, and that the former is measured in bits. 

the mean of the Gaussian distributed Ji in the spin glass 
system (for all i). In order to calculate the concurrence 
for the XY spin glass model, we perform configurational 
averaging over 10^ random realizations. From Fig. 1, we 
notice that the concurrence for the ordered XY model 
vanishes for \i — j\ ^ 3 while in the disordered system, 
it goes to zero for |i — j| ^ 4, implying no significant 
enhancement or deterioration due to randomness. The 
lines in the panels in Fig. 1, connecting the squares for 
the spin glass system and those connecting the circles for 
the ordered system, are the respective exponential fits 
(of Eq.(lO)). To compare entanglement lengths between 
systems with and without disorder, we find that for ex¬ 
ample, with {J)/h = 0.5, ^(^c) = 0.69 in the disordered 
case whereas = 0.50 for the ordered one. The numer¬ 
ical simulations seem to indicate the following: Away 
from criticality, the values of concurrence, for all pairs 
(j, j), are higher in the disordered system as compared to 
the corresponding ordered system, signalling order-from- 
disorder [14, 29, 30] as depicted in Figs. 1(a) and 1(d). 
The roles are reversed as we approach the quantum crit¬ 
ical point (see Fig. 1(b)), except when the factorization 
point is also nearby (see Fig. 1(c)) (cf. [50, 51]). 

Note that although the panels in Fig. 1 are plotted 
for system size = 50, we have checked that increasing 
system size does not change the behavior of the entan¬ 
glements, and hence the entanglement lengths for both 
disordered and ordered systems. 

2. Discord length 

We now show that quantum discord length behaves 
in a qualitatively different way than the entanglement 
length both in ordered and disordered systems. We will 
again compare Di j of the ordered system, governed by 
the Hamiltonian H{J,h), with the quenched averaged 
quantum discord, (Dij), for the XY spin glass system 


with (J) = J/h. To investigate discord length, we di¬ 
vide the entire range of {J)/h (for the spin-glass system, 
which is the same as J/h for the ordered system) into 
three portions, namely 0 < {J)/h < Ac, {J)/h > Ac and 
the neighborhood of A/, with Ac and A/ being the quan¬ 
tum critical and the factorization points of the ordered 
system respectively. Since we investigate the system by 
varying {J)/h, the factorization point lies always in the 
second region and hence the behavior of discord length 
in the second region discussed below are excluding the 
neighborhood of factorization point. The whole discus¬ 
sion will be carried out for 7 = 0.5 and hence the fac¬ 
torization point is A/ = 1.1547. The quantitative feature 
remain unchanged for other values of 7 . 

Case when (J)/h < Ac = 1: In the ordered system, 
quantum discord Dij, decays exponentially as in Eq. (11) 
[26, 27]. In the XY spin glass model, quantum discord 
also decays exponentially, but with a different decay rate. 
As an exemplary case, let us consider {J)/h = 0.5 for 
which the quenched averaged quantum discord of the XY 
spin glass model behaves as 

(A,j) = a + fee”^ ( 12 ) 

with a = 1.4 X 10“^, b = 0.20, = 1.36, while 

a = 4.1 X 10“^, b = 0.18, = 0.56, for the ordered XY 

model, implying = 2.4^£>. Therefore, unlike entan¬ 
glement, we observe significant enhancement of discord 
length in the disordered system as also depicted in Fig. 2. 
The increment of length by introducing randomness in 
the system can therefore only be viewed for quantum cor¬ 
relation measures which are different from entanglement. 
Moreover, we find that in this region, {Di J) > Di j, ex¬ 
hibiting thereby an order-from-disorder phenomenon. 

Case when (J)/h > Ac = 1: Before discussing the 
disordered system, let us first consider the ordered XY 
model. In this antiferromagnetic phase, quantum dis¬ 
cord saturates to a constant value and hence indicates 
long-range order in the system even in the thermody¬ 
namic limit as also predicted in Refs. [26, 27]. For ex¬ 
ample, if one fixes J/h as 1.5, quantum discord of the 
XY model without disorder behaves as in Eq. (11) with 
a = 0.093, b = 0.115 and = 0.80. Note that we chose 
J/h = 1.5 since we are interested in the behavior of dis¬ 
cord length which is far from A/, which is 1.1594 in this 
case. 

In the XY spin glass model, quenched averaged quan¬ 
tum discord again shows long range order. Specifically, 
after an initial decay, it saturates to a constant value. We 
find ^^£ 1 ) = 1.21 > = 0.80 for {J)/h = 1.5. However, 

in this regime, the order from disorder phenomenon is ab¬ 
sent since the value of quantum discord of the disordered 
cases are always lower than that of the corresponding 
ordered ones (see Fig. 2(d)). 

Neighborhood of A/ = ^1 — 7 ^: At the factoriza¬ 
tion point, quantum discord remains constant for all pairs 
of (i, j) for the ordered system and hence goes to 00 , 
for all non-trivial b. In contrast, discord length is finite 
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in the disordered case. Therefore, at the factorization 
point, > C(D>- The decrement of discord length due 
to disorder can also be observed in the vicinity of the fac¬ 
torization point. Hence, enhancement of discord length 
is seen in the entire region of {J)/h except at the neigh¬ 
borhood of the factorization point. 

As a by product, both in the ordered and disordered 
systems, we can prove that quantum discord does not fol¬ 
low any monogamy relation [52] for arbitrary N, in any 
of the three regions. Monogamy of quantum correlation 
quantitatively states that among three or more parties, 
if two of them share high amounts of quantum correla¬ 
tion, then they can have only insignificant amounts of 
quantum correlations with others. Quantitatively, for a 
given N party state a bipartite quantum correla¬ 

tion measure Q is said to be monogamous, with the party 
1 as the “nodal observer”, if 

N 

Q{P1:2...n) > ^ Q{Pli): 

i^2 

where Q{pi: 2 ...n), and Q{pii) are respectively the quan¬ 
tum correlation in the 1 : rest and the 1 : i bipartition. 

Suppose, if possible, that the quantum discord of the 
ground state satisfy monogamy. That would imply 

1 > > ... > N'D{pin). 

The first and the last inequalities are due to the fact 
that each local system is a qubit and T’(pi 2 ) > ’Dipis) > 
'D{pin) respectively. As argued, 'D{pin) can tend to 
a non-zero constant as N —>■ oo, both in the ordered 
and disordered systems. Therefore, in the thermody¬ 
namic limit, N'D{pin) —> oo, giving us a contradiction. 
Hence there exists some N party quantum state for suf¬ 
ficiently N for which quantum discord does not satisfy 
the monogamy relation [53, 54]. Since D{pijs[) can have 
a non-zero value for large N ^ it is easy to see that any 
monogamy-type relations would be violated for quantum 
discord for those states for sufficiently large N . In partic¬ 
ular, a similar argument will imply that also the squares 
of quantum discord also cannot be monogamous for these 
states with sufficiently large N [55]. 

It is also interesting to check the behavior of dis¬ 
cord length with different N. For a fixed {J)/h = 0.5, 
quenched averaged quantum discord, (Dij), with respect 
to the lattice distance r = \i — j\, for different system 
sizes is plotted in Fig. 3 for the disordered system. We 
calculate corresponding to each N. We observe that 
the behavior of quantum discord freezes for > 50 and 
hence we can safely assume that the results obtained for 
A = 50 will mimic those of an infinite spin chain. There¬ 
fore, we take = ^(d)- The scaling of the discord 

length is shown in the inset of Fig. 3 and we find that 
the discord length scales as 



FIG. 3. (Color online.) Scalings of discord length for the one¬ 
dimensional quantum XY spin glass. We plot the quenched 
averaged quantum discord on the vertical axis against lattice 
distance on the horizonal one, for different values of A. The 
logarithms, are with base e. In the main figure, the vertical 
axis is in bits, while the horizontal one is in lattice length 
units. In the inset, the unit of the horizontal and vertical axes 
are respectively logarithms of the total numbers of lattice sites 
and of lattice length. 
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FIG. 4. (Color online.) Entanglement length in random field 
quantum XY spin chain vs. ordered XY chain. All consid¬ 
erations, except for the model considered, and for the fact 
that the different panels are now for different values of J/ {h), 
remains the same as in Fig. 1. 


B. Random field quantum XY spin chain 

We now consider an A-site quantum XY spin chain 
with uniform nearest neighbor exchange interactions, 
Ji = J, but with field strengths, hi, randomly chosen 
from i.i.d., the Gaussian probability distributions with 
mean {h) and unit standard deviation. In the correspond¬ 
ing ordered system, hi assumes a constant value, {h), at 
each site. And again we consider the pure symmetric 
ground state. 

The panels of Figs. 4 and 5 illustrate the features of 
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FIG. 5. (Color online.) All considerations here are the sa: 
as in Fig. 4, except that quantum discord is considered inste 
of concurrence, and that the former is measured in ebits. 


concurrence and quantum discord, respectively, for ci 
ferent choices of (h), against r. We find that the phys 
remain qualitatively unchanged in the random field A 
chain if one compares with the XY spin glasss model. 1 
observe that disorder, in general, does not help in estc 
lishing long range entanglement, while it can significantly 
increase the other quantum correlation length, specifi¬ 
cally discord length, almost everywhere in the parameter 
space except near the factorization point. As observed for 
the XY spin glass model, although discord length, gets 
enhanced due to disorder, the value of quenched averaged 
quantum discord decreases in the region J/{h) > 1. On 
the other hand, the value of entanglement increases in 
this region, showing a complementarity between the two 
types of quantum correlation measures. 


V. ENTANGLEMENT LENGTH AND 
DISCORD LENGTH IN XYZ SPIN GLASS 

In the previous section, we studied quantum correla¬ 
tion length of the quantum XY models. In particular, 
our results showed that even though the disorder-driven 
systems are only minimally benefited in terms of the en¬ 
hancement in concurrence length, a noteworthy endow¬ 
ment occurs in discord length. It is natural to inquire 
whether the findings are generic in one dimensional sys¬ 
tems. Specifically, one can extend the analysis to the 
disordered quantum XYZ spin glass (see Eq.( 3)) in or¬ 
der to find the extent to which the robustness of quantum 
correlations against the distance between interacting sites 
are affected further due to the introduction of additional 
^^-interaction, denoted by A (see Eq. (3)). 

As mentioned earlier, the main difficulty in handling 
a generic one-dimensional system is the absence of an 
analytical approach akin to the XY model. Therefore, 
in order to obtain the ground state of the XYZ spin 



r 
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FIG. 6 . (Golor online.) The ordered and disordered XYZ 
models. Scaling of entanglement as a function of distance for 
the ordered and disordered spin chains with N — 24 and 7 = 
0.5 for {J)/h = (a) 0.5 and (b) 1.5. The up and down triangles 
correspond to the ordered system with A = 0.1 and A = 0.5, 
respectively. The circles and the squares correspond to the 
disordered system for A = 0.1 and A = 0.5, respectively. The 
solid and dotted lines show exponential fits for different cases. 
For the disordered case, the number of random realizations 
taken is 8000. 


chain, we employ the DMRG technique, which is best 
suited for studying the spin chains with open boundary 
conditions in order to achieve high accuracy. However, 
the drawback is that the measurement of the observables 
on the fringes would experience boundary effects. In or¬ 
der to investigate the correlations encapsulated between 
two sites, we consider the central spin on the (A/2)*^ 
site and another site which is positioned at a distance 
r from the (A/2)‘^ site, but is still far from the boundary. 

Let us first concentrate on the entanglement length 
in the ordered as well as disordered systems. Fig. 6 il¬ 
lustrates the behavior of concurrence as a function of 
distance. We notice that in the ordered system, increas¬ 
ing A raises the concurrence length for {J)/h < 1 while 
there is no notable change in the length, when {J)/h > 1 
(see Fig 6). The situation is true also when disorder in 
introduced in the system. However, the disordered sys¬ 
tem fares better than the ordered system in the region 
{J)/h > 1 compared to that of the {J)/h < 1 for higher 
values of A. 

Similar to the XY model, quantum discord behaves 
































quite differently than entanglement in the XYZ model 
with small values of zz- interaction, i.e., A, essentially 
mimicking the results obtained for the XY model. 
Enhancements of discord length are observed both ii“ 
{J)/h < 1 and {J)/h > 1 regions in the presence o 
disorder. For example, for A/h = 0.1, and {J)/h = 0.5 
i{D) = 1-26 for the disordered system while = 0.6- 
for the ordered XYZ spin chain. For higher values o 
{J)/h = 1.5, say, = 1.26 in the case of disorderec 
system and = 1.04 for the ordered one. 

However, in the region (J )/h < 1, the advantage of dis 
cord length obtained in the XY spin glass over the corre 
spending ordered system is faded out with increase of tb 
^^-interaction. It also indicates that there exists an in 
terplay between the randomness in the coupling strengtl 
and the non-random ^^-interaction. In particular, fo 
fixed {J)/h, with the increase of A, the discord lengtl 
decreases in the disordered system. This is illustarted ii 
Table I, where we consider a chain of 24 spins. 



> 

II 

o 

A/h = 0.5 

{J}/h = 0.5 

in = 0.64 
i(D) = 1.26 

= 4.05 

i{D) = 0.86 

{J)/h = 1.5 

in = 1.04 
i(D) = 1.26 

= 0.68 

Ud) = 0.73 


TABLE I. Comparison of discord length for both the orderei 
and disordered systems in quantum XYZ model for differen 
values of {J)/h and A/h with N = 24. 


VI. CONCLUSION 

In summary, we concentrate mainly on the effects of 
disorder in quantum correlation lengths of spin chains. 
The paradigmatic models that we consider are XY 
spin glass model in which coupling strengths are chosen 
randomly, the XY spin chain with random field, and 
the XYZ model with random xx- and yy-couplings. 
Specifically, we compare quantum correlation lengths, 
namely entanglement and quantum discord lengths, of 
the disordered models with the corresponding ordered 
systems. 

We find that entanglement length neither shows sig¬ 
nificant reduction nor increase with the introduction of 
disorder in the system. In sharp contrast, the discord 
length is significantly higher in the XY disordered model 


in comparision to the corresponding ordered ones. The 
features remain unaltered in the XYZ model for low 
values of the z^-interactions. Higher values of the zz- 

:„x_ 1.1 _ 1_ _ _ j.i _i_. _ixu xU - 



r 



FIG. 7. (Color online.) All considerations here are the same 
as in Fig. 6, except that quantum discord is considered instead 
of concurrence, and that the vertical axes are measured in 
bits. 

sites, we prove that quantum discord cannot satisfy a 
large class of monogamy relations for the ground states 
of the ordered and disordered spin models for sufficiently 
large systems. 
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